High-Density Linkage Maps from Japanese Rice japonica Recombinant Inbred Lines Using Genotyping by Random Amplicon Sequencing-Direct (GRAS-Di).

The genetic dissection of agronomically important traits in closely related Japanese rice cultivars is still in its infancy mainly because of the narrow genetic diversity within japonica rice cultivars. In an attempt to unveil potential polymorphism between closely related Japanese rice cultivars, we used a next-generation-sequencing-based genotyping method: genotyping by random amplicon sequencing-direct (GRAS-Di) to develop genetic linkage maps. In this study, four recombinant inbred line (RIL) populations and their parents were used. A final RIL number of 190 for RIL71, 96 for RIL98, 95 for RIL16, and 94 for RIL91 derived from crosses between a common leading Japanese rice cultivar Koshihikari and Yamadanishiki, Taichung 65, Fujisaka 5, and Futaba, respectively, and the parent plants were subjected to GRAS-Di library construction and sequencing. Approximately 438.7 Mbp, 440 Mbp, 403.1 Mbp, and 392 Mbp called bases covering 97.5%, 97.3%, 98.3%, and 96.1%, respectively, of the estimated rice genome sequence at average depth of 1× were generated. Analysis of genotypic data identified 1050, 1285, 1708, and 1704 markers for each of the above RIL populations, respectively. Markers generated by GRAS-Di were organized into linkage maps and compared with those generated by GoldenGate SNP assay of the same RIL populations; the average genetic distance between markers showed a clear decrease in the four RIL populations when we integrated markers of both linkage maps. Genetic studies using these markers successfully localized five QTLs associated with heading date on chromosomes 3, 6, and 7 and which previously were identified as Hd1, Hd2, Hd6, Hd16, and Hd17. Therefore, GRAS-Di technology provided a low cost and efficient genotyping to overcome the narrow genetic diversity in closely related Japanese rice cultivars and enabled us to generate a high density linkage map in this germplasm.

Patterns of genomic differentiation between two Lake Victoria cichlid species, Haplochromis pyrrhocephalus and H. sp. 'macula'.

BACKGROUND: The molecular basis of the incipient stage of speciation is still poorly understood. Cichlid fish species in Lake Victoria are a prime example of recent speciation events and a suitable system to study the adaptation and reproductive isolation of species. RESULTS: Here, we report the pattern of genomic differentiation between two Lake Victoria cichlid species collected in sympatry, Haplochromis pyrrhocephalus and H. sp. 'macula,' based on the pooled genome sequences of 20 individuals of each species. Despite their ecological differences, population genomics analyses demonstrate that the two species are very close to a single panmictic population due to extensive gene flow. However, we identified 21 highly differentiated short genomic regions with fixed nucleotide differences. At least 15 of these regions contained genes with predicted roles in adaptation and reproductive isolation, such as visual adaptation, circadian clock, developmental processes, adaptation to hypoxia, and sexual selection. The nonsynonymous fixed differences in one of these genes, LWS, were reported as substitutions causing shift in absorption spectra of LWS pigments. Fixed differences were found in the promoter regions of four other differentially expressed genes, indicating that these substitutions may alter gene expression levels. CONCLUSIONS: These diverged short genomic regions may have contributed to the differentiation of two ecologically different species. Moreover, the origins of adaptive variants within the differentiated regions predate the geological formation of Lake Victoria; thus Lake Victoria cichlid species diversified via selection on standing genetic variation.

Inferring the demographic history of Japanese cedar, Cryptomeria japonica, using amplicon sequencing.

The evolution of a species depends on multiple forces, such as demography and natural selection. To understand the trajectory and driving forces of evolution of a target species, it is first necessary to uncover that species' population history, such as past and present population sizes, subdivision and gene flow, by using appropriate genetic markers. Cryptomeria japonica is a long-lived monoecious conifer species that is distributed in Japan. There are two main lines (omote-sugi and ura-sugi), which are distinguished by apparent differences in morphological traits that may have contributed to their local adaptation. The evolution of these morphological traits seems to be related to past climatic changes in East Asia, but no precise estimate is available for the divergence time of these two lines and the subsequent population dynamics in this species. Here, we analyzed the nucleotide variations at 120 nuclear genes in 94 individuals by using amplicon sequencing in combination with high-throughput sequencing technologies. Our analysis indicated that the population on Yakushima Island, the southern distribution limit of C. japonica in Japan, diverged from the other populations 0.85 million years ago (MYA). The divergence time of the other populations on mainland Japan was estimated to be 0.32 MYA suggesting that the divergence of omote-sugi and ura-sugi might have occurred before the last glacial maximum. Although we found modest levels of gene flow between the present populations, the long-term isolation and environmental heterogeneity caused by climatic changes might have contributed to the differentiation of the lines and their local adaptation.

Visual adaptation in Lake Victoria cichlid fishes: depth-related variation of color and scotopic opsins in species from sand/mud bottoms.

BACKGROUND: For Lake Victoria cichlid species inhabiting rocky substrates with differing light regimes, it has been proposed that adaptation of the long-wavelength-sensitive (LWS) opsin gene triggered speciation by sensory drive through color signal divergence. The extensive and continuous sand/mud substrates are also species-rich, and a correlation between male nuptial coloration and the absorption of LWS pigments has been reported. However, the factors driving genetic and functional diversity of LWS pigments in sand/mud habitats are still unresolved. RESULTS: To address this issue, nucleotide sequences of eight opsin genes were compared in ten Lake Victoria cichlid species collected from sand/mud bottoms. Among eight opsins, the LWS and rod-opsin (RH1) alleles were diversified and one particular allele was dominant or fixed in each species. Natural selection has acted on and fixed LWS alleles in each species. The functions of LWS and RH1 alleles were measured by absorption of reconstituted A1- and A2-derived visual pigments. The absorption of pigments from RH1 alleles most common in deep water were largely shifted toward red, whereas those of LWS alleles were largely shifted toward blue in both A1 and A2 pigments. In both RH1 and LWS pigments, A2-derived pigments were closer to the dominant light in deep water, suggesting the possibility of the adaptation of A2-derived pigments to depth-dependent light regimes. CONCLUSIONS: The RH1 and LWS sequences may be diversified for adaptation of A2-derived pigments to different light environments in sand/mud substrates. Diversification of the LWS alleles may have originally taken place in riverine environments, with a new mutation occurring subsequently in Lake Victoria.

Cis- and Trans-regulatory Effects on Gene Expression in a Natural Population of Drosophila melanogaster.

Cis- and trans-regulatory mutations are important contributors to transcriptome evolution. Quantifying their relative contributions to intraspecific variation in gene expression is essential for understanding the population genetic processes that underlie evolutionary changes in gene expression. Here, we have examined this issue by quantifying genome-wide, allele-specific expression (ASE) variation using a crossing scheme that produces F1 hybrids between 18 different Drosophila melanogaster strains sampled from the Drosophila Genetic Reference Panel and a reference strain from another population. Head and body samples from F1 adult females were subjected to RNA sequencing and the subsequent ASE quantification. Cis- and trans-regulatory effects on expression variation were estimated from these data. A higher proportion of genes showed significant cis-regulatory variation (∼28%) than those that showed significant trans-regulatory variation (∼9%). The sizes of cis-regulatory effects on expression variation were 1.98 and 1.88 times larger than trans-regulatory effects in heads and bodies, respectively. A generalized linear model analysis revealed that both cis- and trans-regulated expression variation was strongly associated with nonsynonymous nucleotide diversity and tissue specificity. Interestingly, trans-regulated variation showed a negative correlation with local recombination rate. Also, our analysis on proximal transposable element (TE) insertions suggested that they affect transcription levels of ovary-expressed genes more pronouncedly than genes not expressed in the ovary, possibly due to defense mechanisms against TE mobility in the germline. Collectively, our detailed quantification of ASE variations from a natural population has revealed a number of new relationships between genomic factors and the effects of cis- and trans-regulatory factors on expression variation.

A Generalized Linear Model for Decomposing Cis-regulatory, Parent-of-Origin, and Maternal Effects on Allele-Specific Gene Expression.

Joint quantification of genetic and epigenetic effects on gene expression is important for understanding the establishment of complex gene regulation systems in living organisms. In particular, genomic imprinting and maternal effects play important roles in the developmental process of mammals and flowering plants. However, the influence of these effects on gene expression are difficult to quantify because they act simultaneously with cis-regulatory mutations. Here we propose a simple method to decompose cis-regulatory (i.e., allelic genotype), genomic imprinting [i.e., parent-of-origin (PO)], and maternal [i.e., maternal genotype (MG)] effects on allele-specific gene expression using RNA-seq data obtained from reciprocal crosses. We evaluated the efficiency of method using a simulated dataset and applied the method to whole-body Drosophila and mouse trophoblast stem cell (TSC) and liver RNA-seq data. Consistent with previous studies, we found little evidence of PO and MG effects in adult Drosophila samples. In contrast, we identified dozens and hundreds of mouse genes with significant PO and MG effects, respectively. Interestingly, a similar number of genes with significant PO effect were detect in mouse TSCs and livers, whereas more genes with significant MG effect were observed in livers. Further application of this method will clarify how these three effects influence gene expression levels in different tissues and developmental stages, and provide novel insight into the evolution of gene expression regulation.

Whole-genome sequencing reveals small genomic regions of introgression in an introduced crater lake population of threespine stickleback.

Invasive species pose a major threat to biological diversity. Although introduced populations often experience population bottlenecks, some invasive species are thought to be originated from hybridization between multiple populations or species, which can contribute to the maintenance of high genetic diversity. Recent advances in genome sequencing enable us to trace the evolutionary history of invasive species even at whole-genome level and may help to identify the history of past hybridization that may be overlooked by traditional marker-based analysis. Here, we conducted whole-genome sequencing of eight threespine stickleback (Gasterosteus aculeatus) individuals, four from a recently introduced crater lake population and four of the putative source population. We found that both populations have several small genomic regions with high genetic diversity, which resulted from introgression from a closely related species (Gasterosteus nipponicus). The sizes of the regions were too small to be detected with traditional marker-based analysis or even some reduced-representation sequencing methods. Further amplicon sequencing revealed linkage disequilibrium around an introgression site, which suggests the possibility of selective sweep at the introgression site. Thus, interspecies introgression might predate introduction and increase genetic variation in the source population. Whole-genome sequencing of even a small number of individuals can therefore provide higher resolution inference of history of introduced populations.

Factors underlying natural variation in body pigmentation of Drosophila melanogaster.

Molecular mechanisms underlying standing genetic variation of an ecologically relevant trait such as pigmentation trait variation in a model insect, Drosophila melanogaster, are relevant to our understanding of different kinds of intergenomic interactions. In this study, we focused on the association between body pigmentation and stress resistance, and on genotype-by-environment interaction, both of which are likely to contribute to the persistence of phenotypic variation in a natural population. First, we detected a significant association between pigmentation traits in females and starvation resistance (darker strains were weaker) and a weak association between pigmentation and chill coma recovery time (darker strains showed shorter recovery time) among 20 inbred strains from the Drosophila melanogaster Genetic Reference Panel (DGRP), which originated from a natural population in North America. These associations revealed a complex relationship between body pigmentation and physiological traits that may give rise to balanced selective forces acting on the traits under fluctuating environmental conditions. Second, using four of the DGRP strains, a substantial degree of genotype (strain) × environment (rearing temperature) interaction was detected among expression levels of the genes encoding effector enzymes in the melanin biosynthesis pathway. These interactions can potentially reduce the efficiency of purifying selection on the pigmentation traits over a wide range of temperature conditions. Finally, we discuss possible mechanisms that contribute to the maintenance of the standing pigmentation variation in this species.

Complex patterns of cis-regulatory polymorphisms in ebony underlie standing pigmentation variation in Drosophila melanogaster.

Pigmentation traits in adult Drosophila melanogaster were used in this study to investigate how phenotypic variations in continuous ecological traits can be maintained in a natural population. First, pigmentation variation in the adult female was measured at seven different body positions in 20 strains from the Drosophila melanogaster Genetic Reference Panel (DGRP) originating from a natural population in North Carolina. Next, to assess the contributions of cis-regulatory polymorphisms of the genes involved in the melanin biosynthesis pathway, allele-specific expression levels of four genes were quantified by amplicon sequencing using a 454 GS Junior. Among those genes, ebony was significantly associated with pigmentation intensity of the thoracic segment. Detailed sequence analysis of the gene regulatory regions of this gene indicated that many different functional cis-regulatory alleles are segregating in the population and that variations outside the core enhancer element could potentially play important roles in the regulation of gene expression. In addition, a slight enrichment of distantly associated SNP pairs was observed in the ~10 kb cis-regulatory region of ebony, which suggested the presence of interacting elements scattered across the region. In contrast, sequence analysis in the core cis-regulatory region of tan indicated that SNPs within the region are significantly associated with allele-specific expression level of this gene. Collectively, the data suggest that the underlying genetic differences in the cis-regulatory regions that control intraspecific pigmentation variation can be more complex than those of interspecific pigmentation trait differences, where causal genetic changes are typically confined to modular enhancer elements.

The diversity of male nuptial coloration leads to species diversity in Lake Victoria cichlids.

The amazing coloration shown by diverse cichlid fish not only fascinates aquarium keepers, but also receives great attention from biologists interested in speciation because of its recently-revealed role in their adaptive radiation in an African lake. We review the important role of coloration in the speciation and adaptive evolution of Lake Victoria cichlids, which have experienced adaptive radiation during a very short evolutionary period. Mature male cichlids display their colors during mate choice. The color of their skin reflects light, and the reflected light forms a color signal that is received by the visual system of females. The adaptive divergence of visual perceptions shapes and diverges colorations, to match the adapted visual perceptions. The divergence of visual perception and coloration indicates that the divergence of color signals causes reproductive isolation between species, and this process leads to speciation. Differences in color signals among coexisting species act to maintain reproductive isolation by preventing hybridization. Thus, the diversity of coloration has caused speciation and has maintained species diversity in Lake Victoria cichlids.

Correlation between nuptial colors and visual sensitivities tuned by opsins leads to species richness in sympatric Lake Victoria cichlid fishes.

Reproductive isolation that prevents interspecific hybridization between closely related coexisting species maintains sympatric species diversity. One of the reproductive isolations is mate choice based on color signals (breeding color perceived by color vision). This is well known in several animal taxa, yet little is known about its genetic and molecular mechanism. Lake Victoria cichlid fishes are thought to be an example of sympatric species diversity. In the species inhabiting different light environments in rocky shore, speciation by sensory drive through color signals has been proposed by analyses of the long wavelength-sensitive (LWS) opsin gene and the male nuptial coloration. However, the genetic and molecular mechanism of how diversity of sympatric species occurring in the same habitat is maintained remains unknown. To address this issue, we determined nucleotide sequences of eight opsins of six sympatric species collected from a sandy-muddy shore--an ideal model system for studying sympatric species. Among eight opsins, the LWS and RH1 alleles were diversified and one particular allele is dominant or fixed in each species, and we propose that this is due to natural selection. The functions of their LWS alleles were also diversified as shown by absorption measurements of reconstituted visual pigments. To analyze the relationship between nuptial coloration and the absorption of LWS pigments, we systematically evaluated and defined nuptial coloration. We showed that the coloration was species specific with respect to hue and significantly differentiated by the index values of hue (dominant wavelength: λ(d)). The λ(d) value of the male nuptial coloration correlated with the absorption of LWS pigments from all the species, suggesting that reproductive isolation through mate choice using color signals may prevent sympatric interspecific hybridization, thereby maintaining the species diversity in sympatric species in Lake Victoria.

Speciation through sensory drive in cichlid fish.

Theoretically, divergent selection on sensory systems can cause speciation through sensory drive. However, empirical evidence is rare and incomplete. Here we demonstrate sensory drive speciation within island populations of cichlid fish. We identify the ecological and molecular basis of divergent evolution in the cichlid visual system, demonstrate associated divergence in male colouration and female preferences, and show subsequent differentiation at neutral loci, indicating reproductive isolation. Evidence is replicated in several pairs of sympatric populations and species. Variation in the slope of the environmental gradients explains variation in the progress towards speciation: speciation occurs on all but the steepest gradients. This is the most complete demonstration so far of speciation through sensory drive without geographical isolation. Our results also provide a mechanistic explanation for the collapse of cichlid fish species diversity during the anthropogenic eutrophication of Lake Victoria.